


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1984 


Real time prediction of ship hull vibration. 


Taylor, James Wallace 


Massachusetts Institute of Technology 


http://ndl.handle.net/10945/19526 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
. (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist L Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


nN i KNOX appointed — and published -- scholarly author. 
| inp 
, LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


James Wallace Taylor 


REAL TIME PREDICTION OF SHIP HULL VIBRATION. 








JAMES WALLACE TAYLOR 


REAL TIME PREDICTION OF SHEP BULL 


VIBRATION 


Course XIII-A June 1984 








DUDLEY Knox LIBRARY 
na POSTGRADUATE SCHOOL 
ONTEREY, CALIFORNIA 93943-5002 


REAL TIME PREDICTION OF SHIP HULL VIBRATION 
by 
JAMES WALLACE TAYLOR 


B.S., Purdue University 
C1972) 


Submitted to the Department of 
Ocean Engineering 
in Partial Fulfillment of the 
Requirements of the Deqrees of 
OCEAN ENGINEER 
and 
MASTER OF SCIENCE IN NAVAL ARCHITECTURE 
at the 
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
June 1984 
James W. Taylor 1984 


The author hereby grants to M.I.T. and The U.S. Government 


permission to reproduce and to distribute copies of this 
_—-F Lo 


thesis document in whole or Sain — 
« = ‘i 





REAL TIME PREDICTIOCN OF SHIP HULL VIBRATIONS 
by 
JAMES WALLACE TAYLOR 
Submitted to the Department of Ocean Engineering on 
May 11, 1784 in partial fulfillment of the requiements 
for the degrees of Ocean Engineer and Master of Sience 


im Naval Architecture. 


ABSTRACT 


State space models describing the enerqy spectrum of 
the sea, the surface of the sea over an area, anda ship 
hull are developed. These models are used to find a Kalman 
Filter that will estimate the deflection of the ship hull 
from noisy measurements of several points along the hull. 
Various models for predicting the future deflections of the 
hull are tested. 
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Battery alignment is one of the critical factors in 
the control of the weapons systems on board naval warships. 
The spatial relations of the guns, directors, and launchers 
are determined by reference to lines and planes fixed with 
relation to the hull and the determination of these 
relations constitutes the battery alignment. This alignment 
only compensates for the static relations between the 
system elements, and so, when these planes and lines 
distort in a seaway, the alignment is no longer correct. 

Normally, the reference lines used for the aliaqnment 
are the ship centerline for train and the plane of the 
primary element, usually the director, for elevation. When 
the aliqnment is done, the ship is alongside a pier in a 
Knewn condition of laading and in calm seas. At that time 
the relations, vertical, horizontal, and transverse, 
between the directors and the weapons are measured and 
recorded, and the fire control system is adjusted so it 
compensates for static parallax when making the fire 
control solution. Once the ship goes to sea, these 
relations change as fuel is burned, stores used and the 
ship moves in waves. The sea, in addition to causing rigid 
body motions of the ship hull in heave, pitch, surqe, sway, 
roll and yaw, also causes it to flex in the vertical and 


horizontal planes and to twist about its longitudinal axis. 





This flexing distorts the referance lines established in 
the battery aliqnment and introduces errors in the 
targeting solution arrived at by the fire control system. 
These errors are due to dynamic parallax, resulting from 
relative target location due to ship attitude and flexure. 
Studies by the Naval Sea Systems Command, under 
contract to Rockwell International C1], have shown that the 
flexure of the ship hull of monohull vessels appears to be 
significant, relative to riqid body motion. It seems that a 
methed to measure and compensate for these errors is 
mecessary to improve the fire control solution and so to 
increase the fighting effectiveness of the ship. 
Measurement of the flexure could be done directly, for 
example By using a system of sighting instruments, and the 
results fed to the fire control system as they happen. 
Suitable instruments might be Lasers, light beams or strain 
gauges. All of these require the primary equipment, along 
with its support systems, which are subject to failure and 
in addition the system must be aligned and Kept in 
alignment. This results in an increased maintenance load on 
ship’s personel, and the reliability problems inherent in 
mechanical devices. In addition, this method does not 
provide an easy capability to predict what the flexure 
miaqht be at some future time, thus allowing correction of 
the fire control solution at the actual instant of fire. 


An alternative method involves use of a filter ta 





extract deflection information from parameters already 
measured on board the ship. The ship’s gyroscope provides 
ship heading and speed information and in addition it can 
provide heave, pitch, roll, Sway, and yaw information. 
Information regarding wind speed and direction is also 
aVailable. Additional information, if needed, might be 
found from accelerometers placed at the bow, the stern or 
other critical locations. The advantage of using solid 
state devices, instead of directly measuring the 
deflections, is increased reliability and descreased 
maintainace. In addition, the filter formulation lends 
itself to prediction of future distortions, which can be 
fed to the fire control system. 

The objective of this project is the formulation of a 
mathematical model of the ship/sea system that can be used 
with a filter, the formulation of a filter ta estimate the 
deflections, and a predictor to estimate future 
distortions. The equations discribing the ship can be 
written as a set of first order differential equations, 
which lend themselves to a state-space model formulation. 
It is then natural to use a Kalman filter for the 
estimator. This method has already been used to estimate 
the rigid-body motions of a ship for input to an automatic 


landing system [(2,3,4]. 





CHAPTER 1: FORMULATION 


1.1 Literature Overview 


Wave induced ship vibration is the subject of current 


interest to many investigators. The beginnings of the 
majority of the recent work dates from the late 19768’s, 
when the performance of large bulk carriers was 


investigated by Aertssen [5], Aertssen and De Lembre [64], 


Bell and Taylor [£7], and Goocman [8]. These papers were 
mainly concerned with the statistical nature of the 
stresses and bending moments and their predictian, 


Aertssen, and Bell and Taylor doing full scale mesurements 
and Goodman taking a theoretical approach. 

The general trend in the theoretical approaches has 
been to use a modal decomposition, as done for example By 
Eitemepe and FPriee, et. al. £97,16,11,12,13]1, Chen £14], 
Jensen and Pedersen [15], and Kagawa and Fujita [14]. This 
has an advantage when working with reqular waves and for 
seperating out the contributions of the various modes. 

The problem of the exciting hydrodynamic forces has 
been addressed by Salvesen, Tuck, and Faltinsen C1?) for 
the rigid body WSN. using strip theory, providing 
linear expressions for the distribution of the wave forces 
ever the hull. Skj@rdal and Faltinsen [15] have developed a 
linear theory for the hydrodynamics of a springing hull. 


Jensen and Pedersen [15] provide a non-linear formulation. 





Betts, Bishop and Price [19] tie the strip theory approach 
together with structural representations. 

The use of state space techniques for the prediction 
of hull vibration does not seem to have been used by any 
researchers as no references were found by the author 


discribing this approach. 


1.2 Ship Model 


There is no established procedure on how to best model 
the ship and the sea in order to predict flexural 
distortions. One major consideration is the size of the 
models which can make the use of state-space techniques 
very cumbersome. For this reason, for example, it is best 
not to use exact structure and hydrodynamic models, such as 
finite element techniques, and instead to try to reach a 
compromise between model size and modeling accuracy. A 
question arises when trying to combine a structural model, 
such as modal representation, with a random sea. The sea, 
because of dispersion, has a wavelength which is a function 
of wave frequency. Even though the random sea may be easily 
represented by Fourier techniques at a specific point, its 
representation at some other point in space would Be the 
result of the propogation of an infinite number of 
frequency components. How to include these phase shift 
effects in a modal analysis is not clear. In addition, the 


computational load may be large due to the number of modes 





required to accurately represent the quasi-static hull 
distortion, given that the natural frequencies of the hull 
are much higher than the wave frequencies containing the 
major part of the wave energy. 

The approach used here is to approximate the hull asa 
quasistatic non-uniform beam. Using simple beam theory, one 
starts by considering a short section, dx, of the hull ‘(see 
figure 1.139. The governing equation of the beam vibrating 
im a plane of symmetry is found By summing the vertical 


forces acting at point “a® 


mx) dx d/dtCdy/dt) Tox, Comoe YCx, t) + UCx, ti +cV/Gx dx 


fix yc) x + GeV/dx dx C1.1) 


or 


MeazageCay/dt>) = ¢ + VU’ Cit.2) 


Similary, summing the moments about "a*® one obtains 


ia Ax sty = -Vcx,t>? Cleo 


Then by using (1.22) in (1.21) one finds 


m d/dt(dy/dt) = f - M”%’ €1.4) 


-1{6- 





aX 


<> 
a 


r{x,t) ax 


Figure 1.1 
Beam Element 
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From simple, linear beam theory, one has the relation 


Pe ely” * (hes? 


Using (1.24) in (1.235 one obtains 


[eaGieCey/dt) = £ - (Ely’’)’”’ Clie o> 


which is the governing equation for vertical vibrations of 
a Beam. 


If damping is allowed, then one has 


m Oaottary/dt) + &b dy/dt + «Ely’’)’’% = f Clear > 
Once one has the governing equation, the next step is 

to express it in a state-space form. To do this, let xi = y 
and x2 = dy/dt. Then with this coordinate transtoarmation, 


(1.26) becomes 


dxi/dt Se 


nee — CE lxl~~)°" Sm- b x2 /7m + Ff 


The next step is to represent the space derivative. 


Using a finite difference approach, one lets h = cx, the 


-{2- 





distance between points on the hull. Then if there are WN 


points, 


Ee= ONS 1h Cee 


where L is the length of the ship. Expanding the space 


derivative, one has 


eee —— see XA Ot Ue UE! KI’? + Ele x1’*7/ C1.14) 


Now each part can be represented seperately by finite 
differences, with EI, El’, and El’’ as input data. 
At the ends of the ship, the conditions of zero shear 


and zero bending moment lead to 


EI x1’’% = @ 
Ca) 
GET xi°“>’ = g 


Bille a ed eh 


Equations (1.297) with boundary conditions (1.29712 can 
be represented to the order desired and will result in the 


formation of a set of equations, written in matrix form as 


-{3- 





where 


dxl/dt 


Al*X1] + BIxE 


(1.12) 
Y = C) *X] 


Pee lg kl Xt gh Coy ce sg MINI XINZD 
xij = x sub J at point i 
Y = vector of responses at the N points. 


A = a ZN x 2N matrix containing the finite 
difference representation. 


B= distributes the fj’s to the proper state 


—-j{q- 





feo Sea Spectrum 

Ocean waves are a function of both the wind speed and 
duration, among other parameters. One of the commonly used 
wave spectra that attempt to account for these two factors 
is the Bretschneider Spectrum, in which the intensity is 
specified by the siqnificant wave height, H, and the 
duration by the modal wave frequency, Wm. This spectrum has 


the following form. 


SCw) = 1.25 H**¥Z Wm¥e"4d expt—-1.25¢Wm/w) 243 Ci.13) 
- WreS 


The movement of a ship through the water causes it to 
experience excitation not at te wave frequency, w, but at 


the encounter frequency, we, where 


we =w - KU costbeta} Cin 14> 


The energy of the waves remains the same, however, so 


that 








StCwoidw = SCwe ) dwe 
Cla oD 
= S€w> dw dwe 
cwe 
therefore 
Stwe)d = Stw) 
dwe €1.16)3 
du 


where the modulus of dwe/dw is used because Swed is single 
sided. 
Triantafyllou and Athans (€2] have found a rational 


approximation to ¢€1.30> in the following form 


Stwed) = 1.25 H¥*¥2 BCalpha) we £¥4 Clie?) 
Aim C1+Cwe/wod F449) 4£%3 


where 
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alpha = U kim costbeta} Cl Sp) 
9 
BCalpha)d = Nev oor (1.197) 
1+2 alpha 
woCalphad = Wm¢li+alpha) (1.29) 
6.84869 


The corresponding causal system is 


H(s) = So S**¥2 (41221) 
C€14+23¢s/wo)+(s/wo) **¥2) #*¥3 


where 
Sso7— lego H**2 BCal pha? C1.22) 
- Wm 
J = 6.787 
By writing ¢1.38) in standard state space torm (see 


chap 2) one has 


dxe/dt = As*¥xs + Bs¥p 
Ci2o) 


YS Csk*Xs 


== 





where p is white noise of int 


0 i 0 

—“WOX*F2L —-ZJIwo g 

As = 0 G G 
8 8 —“WOX*¥Z 

8 0 G 

8 G 8 

Bs’‘= € @ G G 

Ce = {€{ So 0 G 


Since white noise has a 


frequencies 


ensity q and 


0 0 0 
WOX*¥Z 0 0 
i 8 g Civc24a) 
-2Jwo 0 WOX*Z 
G G i 
G —“WOo*#*2 -—~Z2Jwo 
G G6 pi xwox*3] C1.24b) 
G G g ] €1.24c) 


constant spectrum tor al} 


Sp(w) = q Cl .Z23) 
DI 
and 
Sptw) = Hw) **2 Sw) 
= wo*xxs SO weke*qd WwWoxRrxsd €1.26) 
CWORED + WEK*EGD) H¥S 
one may idenitify the intensity as 
gq = pi wox*s SO C1.27) 
-{38- 





and the driving noise is taken as 


p ~ N(B,1) (1.28) 


with pi wox*3 So included in the model B matrix. 





1.4 Sea Model 

-The representation of the sea must take account of the 
dispersion of the waves in space and still allow for their 
random nature. If the relation between driving point 1, pl, 
and distant point 2, p2, is dependent on a *fixed wave 
number, then the sea at p2 will not be properly related to 
pl. To find a relation that will give the correct relation 
between pi and p2, one starts with Laplace’s equation and 


the free surface condition: 


Qxx + Qzz 6 in the fluid 


CAmeae 
Qtt + g Qz = 6 on the free surface 
Then by assuming that the potential can be written as 
Q¢x,z,t? = UCx,t) expt{Kz3 (1.38) 


and substituting (1.41) into the free surface condition, 


one finds 


Att + g@z = CUtt + QgkU)d expt{kz3 inci? 


If the deep water dispersion relation, K = w**¥2/q, is used, 


€1.92) reduces to 


-—278 -_ 





Utt + we*¥2 U = @ 
or 


Wet eo t = @ 


showing that the free surface condition is satisfied as 
expected. Then substituting ¢€1.41) into the Laplace 


equation, one obtains 


Gece Wixx + UttEt = & Cis) 


It is intresting to note that (1.43) may be obtained 
from the wave equation by using the expression for wave 


phase velocity in deep water. 


and 
Wet = €Cxe2 Wxx = B 
from which 
Wtt -— wk*¥2/(wk#2/Q) **¥2 Wxx = G 
or 
We*¥2 Wtt - g*e*¥2 Wxx = 8 
and since w**2 W = -d/dt(dW/dt> 


Wtttt + g**xZ2 Wxx = G 


Equation (1.943) can now be represented by finite 


-71i- 





difference equations in space and as a set of first order 


equations in time as 


dxv/dt = Av¥xXy + Buys 


Class? 
Ee = Cu*Xyv + Duvkys 
where 
Av’ = Cull,ul2,ul3,ul4,u2l1,...,uN4) 
ujyl = Uy at the j-th point 
uj2 = dUj/dt 
ujs = d/dtCdUj/dt) 
uj4 = Pdttd/dtCdUjy/dt>> 
E = vector of wave elevations at the N points 





1.5 Force Model 

The hydrodynamic exciting forces an a ship hull by sea 

waves have been studied for many years. The approach that 

seems the most fruitful for this project is the strip 
 —_— Of the several forms available (€16,17], the method 
of Gerritsma and BeukKelman [19%] will be used here. 

The hydrodynamic exciting force on the ship may be 
discribed as being the result of the difference between the 
deflection of the ship hull! and the wave onmrofile, z=y-w. 
Then Gerritsma and Beukeliman have proposed that the force, 


FCx,t)d), acting on the ship may be written as 


Ox — eee jMacx) Dz | + NC xd) * Dz + rQBCxdz (1.34) 
Dt Dt Dt 
where 
p> ee dx = a Ud 
dt dx dt at dx 


Then, using z=y-w, and writing the torce as 








Fix, t)) = -Hix,t) + Z2¢x,t? (1.35) 
where 
H(x,t) = Ma(x) D®y + N¢x)-U cMa Dy + rgBixdy (1.36) 
Dt= dx t 
Z¢x,t) = Matx) Dow + N¢x)-U cMa Dw + rgBCxdw Gane. 
Dt= dx Dt 


the force can be seperated into the part due to the hull 


motions and the part due to the waves. The part due to the 


-23- 





hull vibrations may be moved to the left hand side of the 
equation and included in the structural model. The wave 
exciting forces then may be written in a form that may be 
used in the state-space formulation. Then the C-matrix and 
D-matrix of the sea model are modified with equation (1.53) 
to output the force instead of the elevation. The output 
matrix of the spectral model is modified to output the 
first and second time derivatives of the wave elevation in 


addition to the elevation itself. 


1.6 Total model 
The models for the sea spectrum, the sea surface, and 
the ship can be combined into one model by identifying the 


output of one with the input of another in cascade fashion. 


OxX/dt = A¥X + Bp 
(1.38) 
Y = CxX 
where 
A} BixCy B1*¥xDv*¥Cs 
A = 6 AY ByxCs 
A 6 As 
BY =f{C @ G Bs’ J 
C = C Cs 8 6 ] 
X’% = € KI’ Xv’ Xs’ ) 
Deo Ot? 
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CHAPTER 2: STATE SPACE TECHNIQUES 


2.1 State Space Models 


In order to apply Kalman Filtering techniques, one 
must formulate the problem in the form of a state-space 
model. This is a set of first order differential equations 
that discribe the dynamics of the system, its inputs, and 
its outputs. 

The state of a system may be defined at any time as a 
minimum (not neccessarily unique? set of parameters which 
must be specified to discribe that system. A state vector 
is a vector composed of these parameters (state variables). 
The mumber of state variables required is equal to the 
order of the system, n. 

For example, a simple mass-spring-damper system is 


discribed by a second order differential equation. 


m*ed/dt(dz(td/dt) + b¥dzCto/dt + K¥z¢t)d) = FC) G21) 


Then one may define the state variables as _ the 
position, z, and the velocity, dz/dt, of the mass. Next let 
x1 = z and x2 = dz/dt. With this change of variables, (2.1) 


may be written as 


dxi/dt x2 


2? 


nN 


dx2/dt = -K/m*x1 - b/m*xS + £/m 





oe im matrix notation 


dxX/dt = AXX + BxF (2.3) 
where 
6 1 
o—— 
—K/m —-b/m 
BY = f{ 0 1 ] 


mors oe - xd x2 >} 


Alternatively, because the state variables are not 
unique, one may define the state variables as the force in 
the spring, K, and the momuntum of the mass, FP. A 
transformation may be found from one to state vector toa the 


other by using 


K = K¥z 


m*dz/dt 


u 
ll 


Gne may also define an input vector U, containing the 
r inputs to the system and an output vector, Y, containing 
the m outputs. The system equations are then a set of n 


first order differential equations 


dxi/dt 


£1¢X,U,t) 


dxn/dt = #n(X,U,t 





or 


dX/dt = F(X,U,t) 


The output equations relate the output vector, 


the state vector and the input vector 


ym = omex,U,t? 


or 


Y = GCX,U, td 


-?/- 


Ze) 





If the system is linear or linearized, one may write 
the functions fi and qi as linear combinations of the 


states and the inputs. 


n r 
dxi/fdt = ; aij*¥xj + : | bijxuj ce. 5? 
J= j=l 
n r 
y¥K = ) CKj*¥xj + ; dk j*uj Za 7? 
= J= 
feel) nd 
K = (€1,m] 


In matrix notation this becomes 


dxX/dt = A¥X + BU 
C2214) 
Y = CxX + DU 
where A iS an nm xX nm matrix discribing the system 


Gynamics, Bis nx r and discribes the distribution of 
inputs, C is an m xn matrix discribing the distribution of 
states to the outputs, and Dis am x r matrix discribing 
the distribution of inputs to outputs. The form  discribed 


by (2.180) is Known as Standard State-space Form. 


Z.2 Kalman Filter 
In this problem the object is to be able to estimate 
the flexural response of the ship hull at time t and -from 


this to predict what the response will be at time ttdt. To 


-73- 





do this, one common method is to use the Kalman Filter. 
This iS a& method which derives an optimal Cin the least 
squares sensed estimator for the state of a system 
represented by a state-space model. 

Assume one has a mathematical model of the system 


under consideration, 


dx/dt 


Axx + Bw 
(2.11) 


where w is white noise of intensity @, the expected value 
of x at t=6, E[x(6)], is x6 = G6, and the covariance of x at 
t=8, ELx¢@>x’(6d], is X@. Assume further that one has } 


measurements, 2, which are related to x by 


Z= Mex + y Cee) 


where v is measurement noise, modeled as white noise of 
intensity R, uncorrelated with the process noise, w. 

Also assume that the system is observable, meaning 
that the initial condition of any state may be found by 
measuring the output for a finite period of time. 


Then one may propose an observer of the form 


Gy/dt = A*y + K#€O z -— Mky ? Geni o) 


y(@> = EL x(H) J] = 8 


K is Known as the Kalman gain matrix. 


-29- 





One may then define the error in the state estimate as 


e=y-x (2.14) 
and the error covariance as 


P= EC ee” |] C2612) 
It may be shown [23] that the equation which controls 


propagation of the error covariance is 


dP/dt = CA-KM)P + PCA-KM)’% + KRK’ + BQB’ Sr tee 


PCB) 


ACG) =X 


Examining the covariance matrix, P, finds that the 


diagonal is made up of the variances of each of the states, 


P= EL ee’ J = el [Mel vez ... en ] 
eZ 
en 
elx*2 elxez : ; A 
e2xel e2**¥Z2 5 ; é 
is . : ; C2207) 
. : enk*ez2 


So if one takes he Trace of P 


-—34a- 





Tr¢P) = e| **2 (2.18) 


one has the sum of the squares of the errors in each of the 


states. Athans [(€23] shows, using the matrix minimum 
prinicipal, that the K that minimizes the Trace of P, ie; 
the sum of the squares of the error, is 

K = FM’R G2 li) 
If, as in the present case, none of the matrices A,B,M,Q 


and R are functions of time, then dP/dt — 8 in the steady 


state solution, and one has 


6 = AP + PA’ + BQB’ —- PM’R MP (2.268) 


Summarizing, then, if one has a linear system 


described by 


dx/dt = A*x + B*w 
(2.21) 
p = Cx 


one may estimate, with minimum square error, the state x by 


using y which is described by the equation 


h 
th 
Naw? 


Ge 7Gee—eney + KC 2—-My) CZ. 


IN 
a8) 
“er 


K = PM’R CZ. 


-3i1- 





with P given by Cin the steady state case? 
S@ = AP + PA’ + BRB’ —- PM’R MP (2.24) 


Software pacKaqges exist, for example CTRLC £26], which will 
solve 2.24 for the steady state values in P and so for the 
values in K. 

Once the values of the state at time t6 are Known, 
either from direct measurement or from estimates, the next 
task is to predict the the states at some future time, 
tf=tG+dt. To do this, one may solve 2.271 subject to 


xCt@>=xB as 


xCt>? = explACt-tGd sxe + exp CAct—s) Bw sdds C2 aa 


tG 


Now, suppose that at time t¥, one wants to predict 


xCtB+dt>. Then 


t 
xCtf> = explAtCdt>>*xa + expfACtf—-s) }Butsdds (2.26) 
t@ 
but té<s<=tf is some future time and ws) is unknown 


because w is random. In fact, since w is white noise and so 


completely random, it is completely unknown and the est 





estimate of it is its average value or expected value, 


EL w(t) J] = @ 


Therefore, the best estimate of x(tf) is 


xpCtf) = exptAdti#*r (ta) 


where y has been used in the initial condition 


is obtained from the Kalman filter. 


eee) 


26 2O0 


because it 





2.3 Approximations 


The question of the order of the finite difference 
representation will be addressed next. The accurarcy of the 
distortion results are dependent on that order. For just 
one point of the ship model, say Jj, equation 1.27 en be 


written as Cwith j understood in the coefficients? 


Oxji/dt = xj2 
(2.2973 
Geperat = —E] xji”“?“* = El’ xyi’’’% = 
mm nm 
Bi x) eae Xie + Ff 
m Tm Tn 


Since the highest derivative in x is the controlling factor 
in the discretation error, the case of a uniform ship, with 
El’ and EI’’ both equal to zero will be investigated. The 
method used will be that used by Smith [23]. 


At point j, €2.30) becomes 


eet = Axx? ~~ + DEX + F (2.30) 


where the matrix A Is 


i g 
oe 
Se ti 6 


ana © is 


6 g 
oo 
8 -b/m 





and F is 


#/m 








Now write the tourth derivative of X as 


pee = AC ItZ) H4XC tid + 6X03 = 4AXC§~-i1)d + ACj—2) 


h*x*4 


Being €2.32) in €2.31), one finds 


dx¢(jo/dt = A XCjt2) -dAXCJ+1) + 6A XCJ) 
h**4 


- 4AXCj-1) + AKCj—-Z2 
h**4 


+ D XJ) 


and expanding the X¢€jti)d by Taylor series 


XCgtid =. + Rex” + {1 WHxex2EXS?" + 


HhE¥SEXS “SZ + 
Zi 


we 


I" 


Substituting (2.34) into (2.33), Keeping terms to 


h*x*8 and collecting terms leads to 


agxcjo)7dt = A¥X’’*’% + DEX + F + 


Asx ehex2 ) +¢ © 1 OR eh exG D 
Bo 


( 1% 
é 


(2.31) 


order 


(2.34) 





where the term in the first bracket can be interpeted as 
the order h**Z error and that in the second bracket as the 
order h**4 error. 

In order to see the effect on a particular solution to 
(2.36), consider the solutions of the natural modes of a 
uniform beam (£22). The space dependent part is composed of 
sinusoids and hyperbolic functions in the variable Kx. 

wx) = sin€Kx)? + sinh(Kx) + s¥¢C cos(kKx) + caoshtKx) ) €2.35) 
where $s is a constant dependent on the mode.If L is the 
length then 
KE ~ Cj+ dF 1 
defines the approximate values of kK. Taking the eight 
dirivative of ¢€2.35) brings out a K*¥*3 term, so 


W = K*¥*¥S¥w 


and (2.31) can be written as 


Grecit —eemn 97" + D*X + F + E (2.345) 
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where 


ee eh ead 


2 


&@ 


AXxK*¥*¥S¥X¥h¥x4 


1 


3G 


Then since 


Kx¥*¥GeX = X°°~"" 


one may write 


GX/dt = AxX7“%" €C 1 + 1 #HRR4EKE4 D + DEX + F 2.57) 
88 


Then the relation between the point spacing and the error 


may be found from the second term in the brackets 
= = 1 #h¥¥4RKEHS 
85 


I# the maximum desired error is eo, then 


hx*¥4 < 88*¥ e0 
Kx*4 
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or 


he S229 Fe eo C2.c5> 
Cite Pi 


By similar means the relation for order h-squared can 


be found to be 


h < 2.45*L eo (2.39) 
Cj+2)Pi 


Table 2.1 shows h vs. eo for beth orders for the 
springing mode, and table 2.2 shows h vs. j for a fixed 


eo=G@.125. Length is 188 meters. 


Table 2.1 


h vs. eo Cj=1) 
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Table 2.2 


ho vs. j Ceo = 6.125) 





The most obvious fact to be gained from Table 2.1 is 
that for a given accuracy, the number af points and hence 
the number of states is much less with the order he¥*4 
approximation. Each added point requires six additional 
states in the total model, two for the ship model and four 
for the sea model. However, only five points are required 
to approximate the fourth direvative to order h*¥*2 but 
eight points are required for order h*¥¥*¥4. These numbers set 
the minimum number of points required for any accuracy. 

As the number of nodes in the vibration increases, the 
number of points for a qiven accuracy increases rapidly, as 
seen in Table 2.2. Aqain, the order h*¥*4 approximation is 
the less restrictive of the two. It only requires & points 
to achieve 12.5% accuracy with order h**4, while 1? points 
are required for the order h*¥*2 approximation. 

This means that while it is possible to use an order 
h*¥*¥2 approximation, an actual inplementation will require 


both a higher order and a greater number of points. 





CHAPTER 3: NUMERICAL RESULTS 


emi otability of Models 

The mathematical models of the sea spectrum and the 
ship structure are stable as modeled. However the force 
model is not. The fourth order equation discribing the sea 
surface has four solutions, two of which are exponentials, 
and it is these solutions that can cause instability. The 


equation is 


Utttt + Qqe*2 Uxx = @ Kot 
which has solutions in the form 
Ucx,td) ~ € exptCiwt), expC-iwt), exptwt), exp C-wt) 3 
CSae? 
The exponential roots, in the analytical problem, are taken 


care of by the initial conditions. When using a computer, 


however, even if the initial conditions are set equal to 
zero subsequent round-off errors introduce fictitious 
non-zero initial conditions, which excite growing 


exponentials. 


To see this, consider the simple example 


mte > a4"42 Y = Y (3.3) 
¥(6) = 6 


YCinFfinity?d G 
which has the analytic solution 


Y¥Ct) = A expCat) + B expat) €3.4) 
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By approximating this by finite differences, one obtains 


Ytt ~ VYCitld -2RYCI)D + YCIi-1) = ax*¥*¥2 YCi) CS. 
a=. ° 
or 
YC1+t1) =-€2 + Cahd*#*2) YCI)D + YCi-1) = &G 
To solve this finite difference equation, take 
YOn) = rx¥n 
then one obtains 


rm¥*¥2 —€2 + Cahd##2) r+ 1 = a 


and then for h <¢ {7a 


Rew oi + ah Fee ah**2) 
and 
See ES ah + OCah*xx2>) 
Then taking a linear combination of the solutions 


YOnd) =A C1 + ahd*¥en + B C1 - ahd ¥¥n (3.6) 


it can be seen that no matter how small A is, if it is not 
exactly equal to zero, the part of the solution involving 
the growing exponenitial will eventually dominate the 
solution, swamping the correct solution. If the objective 
were only to obtain the solution to the equation, one might 
try integrating backwards in time. Here though, the 
objective is to produce the solution in real time, so the 
equation must be integrated forward in time. 

Thie¢ is mot as severe a problem as it might seem at first 


look. One of the properties of the Kalman Filter is the 
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stabilization of the filter model. The gain that results 


will result in a closed loop observer matrix, 


Ac] =A - K*¥C 


that has all its eigenvalues in the left-half plane. The 
only place the problem arises is in simulation of the 
system. By defining a reconstruction error, e=x-xhat, and 
using the expressions for x and xhat, it can be seen that 


the error must obey the differential equation 


geen? = Ea-KCI ¢< {> (3.7) 
dt 
This means that the reconstruction error qoes to zero as t 
becomes larqe only if the closed loop observer matrix is 
asymptotically stable. 

Simulation of the system is necessary to test the 
observer as the model is varied. The method used is to sum 
a number of sinusoids, each of which has an amplitude 
proportional to the energy in the sea spectra at that wave 
frequency and a phase that is random with respect to the 
phases of the other sinusoids. 

Assume that we are given a sea spectrum, Sw), and we 
wish to simulate the sea elevation at a number of points. 
Since the instantanous sea elevation at a point x and 
frequency w can be represented as the square root of the 


energy at that frequency ‘(multiplied by 2 since the sea 
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spectra is only defined from @ to infinity), 


Vex, wd) =VZ SCx,w) 


then the elevation at x, including all frequencies should 


be 


Cu 
aa | 
VUCx) = Ret |\J2 Stw? dw! expt iptw) 33 eee 


8 


where p is the random phase. Then including the time factor 


one has 
CO 
Vix) = Ret \V2 Stw) dw expt itwt + ptw) 33 Coa) 
a 


For a computer simulation, equation ¢€3.%7) must be 


discretized to become 


Nw 
VOx,t) = {V2 Stwj) Awj' cost wit + p(wj) 3 (3.16) 
ao 
where here p is a random number uniformly distributed 


between 6 and two pi. 


The simulated sea surface is then used to drive the 


Ao 





ship model to form a simulation model, the output of which 
is used as a source of the measurements supplied to the 
observer model. This model was run and the results stored 
in a computer file for later use. Figure 3.1 shows a 


representative run of the simulation. 


3.2 Estimation and Prediction with Correct Model 

The Kalman filter for the total model is required to 
construct optimal estimates of the states of the system 
that can be used to predict the future output of the 
system. This includes those states that cannot be measured 
directly, such as those of the sea. To do this it must have 
a source of information about the system. 

The filter gains are affected not only by the model 
used but also by the noise in the measurements themselves. 
This noise results from a number of sources. For example, 
from measuring instrument noise, slamming of the hull, 
local vibration effects and vibrations excited by external 
sources. The filter then must reject the hiqher frequency 
noise to obtain the needed information about the states. 


The measurements are found from 


Zoey tae KON B,D Csr) 


where R is the measurement noise covariance matrix, defined 


as the expected value of the noise process 
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R= EC vv’ 3 (3.12) 


Then what is needed are estimates of the measurement noises 
that would be encountered in an actual ship hull. Using 
studies of ship vibration, the noise covariance . for 


displacement was taken as 


rCi,id = 3.6 E -4 meters squared 


at amidships. According to Lewis [27], the amplitude 
amidships will be equal to about 6.25 the value at the end 
fer machinery itnduced vibration. For velocity of 


displacement a covariance of 


rCi,i>? = 6.6 E -1 meters squared / second squared 


at amidships was used. All off diagonal values were assumed 
to be equal to zero, implying that the measurement points 
are isolated from each other. 

The filter was calculated using the models for the sea 
spectra, forces and ship in cascade. This allowed an 
investigation to determine the required measurements to 
insure the observability of the system. In addition, the 
filter pole locations were determined. For the Great Lakes 
Ore Carrier Stewart J. Cort, described in Appendix A, the 
poles were all located at a frequency of less than 1.85 
radians per second. The pole locations are shown in Table 


3.1 for the output matrix 
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16666 6 6 8B 

C= SBiesesbses se @ CS) 
661666 6 8B 
666164 86 6 Bb 

This output matrix provides measurements OF the 


displacement. The resulting displacement estimation errors 
are given in Table 3.2 and Figures 3.2 and 3.3 show the 
estimated value compared with the simulated value. The 


noisy siqnal is shown tn fiqure 3.4. 
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Prediction of the future values of the 


moee 3.1 
Kalman Filter poles 
Se, 


—~-§.24E-4 + i 2.87 
-7.S@E-3 + i 1.488 

-9.836E-9 + i 7.38E-1 
-6.67E-3 + i 6.46E-1 
-1.5@E-2 + i 4.75E-i1 
—-§.39E-1 + i §6.58E-1 
-3.498E-1 + i 6.22E-1 
-~§.32E-1 + i S.SiE-!1 
—-$./SE-1 + 1 1.33E-2 
-7.37E-1i + i 4.92E-4 
-6.46E-1i1 + 1 4.58E-3 
-5.95E-3 + i 4.86E-2 
—-§.8HE-4 + | 4.&SE-2 

Table 3.2 
RMS Estimation errors 

(meters) 

 .36E-2 

9 .36E-2 

Loo — 1 

Y .s6E-2: 


states, 


qQiven 


the current measured or estimated values can be dane by the 


equation 


let 4 


gets large. 


the future the prediction 


is stable, 


Gxp 
dt 





xp(@) 


A * xp 


xhat (a) 


the predicted value will 


The abjective then 


prediction time step, 


is ac 


the actual 


is to determine how far 


ceptable. If dt 


prediction Becomes 


ic 


(3.14> 


go to zera as time 


on 


the 


xp¢t) = expt A *dt 3} * xhatctd (3.15) 


which involves matrix-vector multiplication only. 

For this system, a number of prediction models were 
tried, but none gave usable predictions. The first model 
tried was to use the unstabilized system, but as expected, 
the output immediately started growing. 

Next, the system model was stabilized by using a very 
high measurement noise to calculate the feedback gain. For 
this model, the predicted deflections showed an initial 
change in the same direction as the true value But rapidly 
decreased in amplitude to about one tenth of its amplitude. 
The phase of the prediction was very close ta that of the 
true deflection. 

Another method tried was to stabilize the madel of the 
sea surface only. This resulted in a stable total mode! But 
no improvement was seen in the prediction. Here the 
predicted value rapidly diverged from the true value after 
about four time steps and the phase showed no apparent 
relation to that of the true phase. 

Finally, a prediction was tried by using only the 
model of the ship. This was even less successful as the 
prediction immediately diverged from the true, often 
increasing when the true value was decreasing and 
vice-versa. 

Though prediction was not acheived by any of these 


models, this does not mean that prediction can not be done. 
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It iS one of the underlying ideas of the state space 
technique that all necessary intarmation about a system is 
included in the state of the system at any time. Prediction 
models that bare investigation are: (1) stabilize the sea 
surface model by treating it as a regulator problem = and 
selectively choosing the elements in the control weighting 
matrix to limit the wave amplitude, €2) try the same method 
on the total model, 3) a useable, if short prediction may 
result from the second method tried if the time step is 
reduced, ¢€4) use the first method tried but scale the 


predicted output. 





CHAPTER 4: DISCUSSION 


The problem addressed in this thesis is the estimation 
of ship hull vibrations from a minimum number of measured 
parameters. These estimations are then used to predict the 
future deflections of the ship hull. 

The first step was to develop a mathematical model of 
the ship-sea system in state-space form. This involved 
three seperate models. The first, which models the sea 
pewer spectrum, was taken from the literature. This mode} 
is based on the ITTC ocean spectrum that uses as parameters 
the modal wave frequency and the siqnificant wave height. 
In this spectra, the modal frequency can be taken as 
4.97/Tave, where Tave is the average period measured at the 
calm water level. The significant wave height may be taken 


#rom table 4.1 if no measurements are available. 


Table 4.i 


wind speed (Knots) sig. wave height (m) 
29 Sei 
36 Jel 
48 Ss. 
28 11.8 
68 14.6 


The model to include the dispersion relation for sez 
WAVES Was developed from a velocity potential 
representation of the sea, resulting in an unstadic model. 


The instability was due to the increasing expanential 
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solution to the resulting fourth order in time differential 
equation. 

The third model models the ship hull itself. This was 
derived by assuming that the ship could be represented as a 
simple beam. Then a composite, continuous time, discrete 
space, finite difference representation of the simple beam 
equation was used for the final model. The model was 
developed using a second order finite difference 
approximation but for an actual applicatian, a fourth order 
approximation should improve the accuracy to a significant 
degree. 

Gerritsma and Beukelman’s relation relating the wave 
height to force on the ship was used to couple the sea 
model with the hull model. Since in a state-space 
formulation the time derivatives of the wave height and the 
Geflection are available in the state vector, implementing 
the force model simply required additional terms in the 
output matrices of the sea spectra and sea surface models, 
and added terms in the hull system matrix. 

Using these models, CTRLC (€26], a software package for 
control system design, was used on a Digital VAX-11/7S2 at 
the Jaint Computer Facility to do the simulations, develap 
the filter gains and evaluate the sensitivity to errors in 
the model. The estimation of the deflections of the hull 
could be done well with the state space model developed. 


The prediction of the deflection to some time in the future 
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was not successful using the models developed here, but 


with improved models it could be possible. 
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APPENDIX 
BREDICTIGN RESULTS 

The several methods of prediction attempted gave less 
than satisfactory results when compared to the simulated 
Seal . Fiqure A.1 shows the prediction that resulted trom 
usinq the system model stabilized By using a very large 
measurement noise. This resulted in a predicted value that 
was much larger than the simulated value. The same general 
pattern resulted when the measurement noise was varied trom 
& low equal ta the expected measurement naise to as high as 
could be computed. Nor did changing the starting point of 
the prediction with relation to the phase ot the simulated 
signal improve the prediction. 


The next tiqures show the prediction resulting from a 


mode} constructed by stabilizing the sea-surface mode} 
berore forming the total model. In Figure A.2Z, the 
predicted signal is very similar to the prediction in 


Figure A.1. If just the first few time steps are plotted, 
it can be seen that the predicted signal does not track the 
Simulated signal even for a short period of time, Fiqure 
eas 

Finally, Fiqure A.4 shows the results of using only 
the ship model to do the prediction. As can be seen, this 
prediction has no apparent relation to the simulated 


signal. 
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